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ABSTRACT 



Context. Many efforts are being made to characterize extrasolar planetary systems and unveil the fundamental mechanisms of planet 
formation. An important aspect of the problem, which remains largely unknown, is to understand how the planet formation process 
depends on the mass of the parent star. In particular, as most planets discovered to date orbit a solar-mass primary, little is known 
about planet formation around more massive stars. 

Aims. To investigate this point, we present first results from a radial velocity planet search around red giants in the clump of 
' (— I ' ' intermediate-age open clusters. We choose clusters harbouring red giants with masses between 1.5 and 4 Mq, using the well-known 

cluster parameters to accurately determine the stellar masses. We are therefore exploring a poorly-known domain of primary masses, 
I . which will bring new insights into the properties of extrasolar planetary systems. 

O ' Methods. We are following a sample of about 115 red giants with the Coralie and HARPS spectrographs to obtain high-precision 

' radial velocity (RV) measurements and detect giant planets around these stars. We use bisector and activity index diagnostics to dis- 

C/3 \ tinguish between planetary-induced RV variations and stellar photospheric jitter. 

. Results. We present the discoveries of a giant planet and a brown dwarf in the open clusters NGC 2423 and NGC 4349, orbiting 

the 2.4 Mo-star NGC2423 No3 (TYC 5409-2156-1) and the 3.9 MQ-star NGC4349 Nol27 (TYC 8975-2606-1). These low-mass 
companions have orbital periods of 714 and 678 days and minimum masses of 10.6 and 19.8 Mjup, respectively. Combined with the 
. ^ other known planetary systems, these detections indicate that the frequency of massive planets is higher around intermediate-mass 

, stars, and therefore probably scales with the mass of the protoplanetary disk. 

'NT . 

. Key words, stars: individual: NGC2423 No3 - stars: individual: NGC4349 Nol27 - stars: planetary systems - galaxy: open clusters 

and associations: individual: NGC 2423 - galaxy: open clusters and associations: individual: NGC 4349 - techniques: radial velocities 

(N 

' 1. Introduction quently around solar - mass stars than around low-mass stars 

(lLaughlinetal.ll2004l:llda&Linll2005i) . This is also supported 
O More than 200 extrasolar planets have been discovered over the observational data showing that giant planets a re much rarer 
L; past decade, revealing a completely unexpected diversity in the ^^ound M dwarfs than around solar-type stars (iBonfils et al.l 
properties of planetary systems dUdry et al. || 2007D . Models of 2005; Endl et al. 2006; Butle r et al. 2006). 
planet formation had to be largely revised to explain a num- 
ber of characteristics derived from the observations such as the I" this context, it is of great interest to investigate the fre- 
planet-metallicity correlation, orbital migration, and the occur- q^ency of giant planets around stars significantly more massive 
rence of high eccentricities. Increasingly complex models in- *an the Sun. Unfortunately, it is impossible to apply usual high- 
cluding Monte C arlo simula tions are being develope d by dif- precision Doppler techniques to stars with spectral types ear- 
ferent groups (e.g. Ida & Lin ' 2004; AIibertetalJ|2od5|) with the li^r than late-F due to the increase in rotational velocities and 
goal of comparing theoretical predictions to the observed distri- *e low number of spectral lines in these stars. A few attempts 
bution of extrasolar planets have been ma de, however, to adapt the technique to A-F stars 

Besides metallicity, protoplanetary disk masses and surface dGalland et al^ . M, unveiling a brow n dwarf candidate around 
densities are other important factors suspected to strongly in- ^ 1.7 M© primary d Galland et al. || 2006|) . 

fluence the formation of giant planets. They are thought to Another promising approach is to study intermediate-mass 
be dependent on stellar mass, in the sense that more massive stars in a more evolved stage, i.e. in the red giant phase, where 
stars wi U have more ma ssive disks and higher surface densi- rotational velocities are low and many spectral fines are avail- 
ties (see llda & LinI 12005]) . altfiough this point needs confirma- able for Doppler shift measurements. However, great care has 
tion. The exact impact of this on planet formation is presendy to be taken when interpreting RV variations in red giants due 
poorly known theoretically but the favoured core-accretion see- to the presence of intrinsic stellar jitter Precise RV measure- 
nario predicts that massive giant planets should form more fre- ments have afieady been reported for many red giants in the so- 

lar neighbourhood, and a few planet candidates have been dis- 

Send offprint requests to: C. Lovis covered around such stars ( Frink et al. 2002'; Hatzes et al.l l2003l 

* Based on observations made with the ESO 3.6m-teIescope at La 2005, 2006; Setiawan et al. 2003, 2005; Sato et al. 2003). These 
Silla Observatory under program IDs 075.C-0140, 076.C-0429, 077.C- detections demonstrate that in these cases steUar jitter does not 
0088 and 078.C-0133. hide the RV signal of giant planets. It appears that a significant 
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fraction of red giants are indeed suitable for planet s earch, pro- 
vided they are chosen not too evolved or too cool dSato et alj 
I2005h . 

In the context of planet formation around intermediate-mass 
stars, it is crucial to have a precise knowledge of primary masses. 
This is unfortunately very difficult for field red giants such as 
those mentioned above, due to the degeneracy of evolutionary 
tracks in the HR diagram. Error bars on the stellar mass are often 
of the order +2Mq, making it impossible to determine whether 
the star under consideration has a mass around 1 Mq or signifi- 
cantly higher. 

To overcome this difficulty we started a high-precision radial 
velocity survey of red giant stars in a number of intermediate- 
age open clusters. Knowledge of cluster parameters allows us 
to accurately determine red giant masses, which we choose be- 
tween 1.5 and 4 Mq. We obt ained first measurem ents in 2003 
with the Coralie spectrograph (Ouelo z et alj|200 01 at the Swiss- 
1.2m Euler telescope at La Silla Observatory, Chile. However, 
the faintness of our targets prevented us from observing all 
of the clusters with Coralie and led us to apply for observing 
time on a more po werful instrument, the HARPS spectrograph 
dMavor et al.ll2003h on the ESO-3.6m telescope, also at La Silla. 
This allowed us to reach a sufficient RV precision on the faintest 
red giants {V = 12) in our sample. 

Ve ry recently, following a strategy similar to ours. lSato et alj 
(l2007h published the discovery of a giant planet orbiting a 
2.7 Mq red giant in the Hyades. This illustrates the potential 
of open clusters as laboratories to understand planet formation. 
In this paper we present the discoveries of a massive planet and 
a brown dwarf around two intermediate-mass stars in our clus- 
ters, NGC2423 No3 (TYC 5409-2156-1) and NGC4349 No 127 
(TYC 8975-2606-1). Sect. 2 describes the sample of stars we are 
following in our RV survey. The basic properties of both planet- 
host stars are discussed in Sect. 3, whereas Sect. 4 presents the 
RV measurements and orbital solutions for both companions. In 
Sect. 5 we show that the RV variations are not due to spurious 
stellar effects. Finally, we discuss the implications of these dis- 
coveries in the broader context of planet formation in Sect. 6. 

2. Description of the survey 

The clusters in our survey have been selected based on a num- 
ber of criteria. First of all, to have giant masses between 1 .5 and 
4 Mq, they have to be of intermediate age (0.2-2 Gyr). They also 
have to be observable with Coralie and HARPS, so we set a dec- 
lination limit 6 < +15°. Moreover, to reach a sufficient photon- 
limited RV precision (<10 m s"'), their giants must be brighter 
than V = 10 for CoraHe and V = 13 for HARPS. Finally, we 
only took into account clusters having more than 3 giants known 
to be non-binary cluster members. 

Table [T] gives the list of all selected clusters with their 
main properties. Cluster ages and data o n their red giants 
have been taken from the W EBDA database ( MermilUodll 19951; 
iMermilliod & PaunzenI 12003) . The mass of the red giants in 
these clusters was derived from their ages u sing the Padova stel - 
lar evolution models at solar metallicity dGirardi et alj l2000l) . 
Error bars on the masses have been estimated taking into ac- 
count the uncertainties on the cluster ages and metallicities. The 
assumption of solar metallicity is not likely to cause large er- 
rors on the derived masses since the metallicity distribution of 
open clusters in the solar neighbourhood shows a peak at solar 
metallicity and a small scatter of ~0.1 dex (Twarog et al. 1997). 
The determination of precise cluster metallicities is a challeng- 
ing task and the values found in the literature for a given clus- 



ter are sometimes quite different. This is especially problematic 
for planet searches since we would like to distinguish between 
two effects, the well-known planet-metallicity correlation and 
the influence of stellar mass on planet formation. For this rea- 
son we are planning to obtain high signal-to-noise spectra of the 
giants in our sample in order to perform a high-precision spec- 
troscopic analysis and de termine met allicities in the same way 
as for nearby field stars dSantos et ai.i ,2004; .Valenti & Fischen 
120051: Ida Silvaet al.ll2006l) . 

3. Parent star characteristics 

3.1. NGC2423No3 

Our sample contains 6 giants belonging to NGC 2423 (distance 
d - 766 pc). With an age of ~750 Myr, these stars have an esti- 
mated mass of 2.4 + O.2M0. The metallicicty g iven in WEBDA 
for this cluster is [Fe/H] = 0. 14 + 0.09 (from Twarog et aljl997l) . 
We verified that all 6 stars are indeed cluster members by check- 
ing their position in the HR diagram and their radial velocities. 
Fig.[I]shows the HR diagram of NGC 2423 as given by WEBDA 
with the position of our 6 survey stars. All are situated in or 
close to the red clump. NGC2423 No3 (a = 07h37m09s, 6 = - 
13°54'24") has visual m agnitude V - 9.45 and colour index 
B-V - 1 .21 according to 'Hassan' ('1976'). It appears to be some- 
what redder and brighter than the mean clump position, indicat- 
ing that it might be slightly more evolved. However, the small 
number of red giants in the cluster and the uncertainties in the 
photometry make it difficult to draw a clear conclusion on its 
precise evolutionary stage. Regarding radial velocities, we com- 
puted from our Coralie measurements a mean radial velocity for 
the red clump of 18.67 km s"', with a dispersion of 0.34 km s"'. 
The mean RV of NGC2423 No3 is 18.32 km s^', leading to a 
high probability that it is indeed a cluster member. 

We started the monitoring of the NGC 2423 giants with 
Coralie in 2003 and accumulated ~10 RV measurements per star, 
spread over more than 2 years. The average visual magnitude 
of these stars, V - 10.5, makes them difficult targets for pre- 
cision RV measurements with this instrument. Nevertheless we 
could obtain a typical photon-limited RV precision of -15 m s 
sufficient to detect giant planets. Fig. |2] shows the histogram 
of the RV dispersion for the NGC 2423 giants after correc- 
tion of the instrumental contribution. Four of them exhibit 
a RV scatter below 30 m s"', while NGC2423 No3 stands 
clearly out with its RV dispersion of 72 m s '. Finally, the last 
star (NGC2423 No43) shows a large-amplitude, long-term drift 
(ctrv = 54 2 m s~') due to a stellar- m ass companion recently dis- 
covered bv lMermilUod et al.l (12007 '). 

These results allow us to compute an estimate of the typical 
RV jitter affecting the NGC 2423 giants. Quadratically subtract- 
ing the instrumental contribution, we obtain a value of ~20 m s"' 
for the jitter, based on the 4 stable stars in this cluster. It therefore 
appears that the jitter level will not prevent us from detecting the 
RV signal of giant planets orbiting these stars. 

3.2. NGC4349No127 

NGC 4349 is a distant open cluster (d = 2200 pc) with an age of 
only 200 Myr. The be st available estirn ate for the metallicity is 
[Fe/H] = -0.12 ± 0.04 (iPiatti et al.lll995l) . This is one of the most 
interesting clusters in our sample since its giant stars have an es- 
timated mass of 3.9 + 0.3Mq. We have included 7 of them in our 
survey, whose positions in the HR diagram are shown in Fig. [3] 
NGC4349 Nol27 (a = 12h24m35s, 6 = -61°49'12") has visual 
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Fig. 1. Colour-magnitude diagram for NGC 2423 taken from 
WEBDA. Red giants in our survey are shown as triangles, 
NGC2423 No3 as a square. 



magnitu de V = 10. 88 and colour index B -V - 1 .46 according 
to Loh mannid 19611) . Again, it seems to be slightly more evolved 
than the mean clump position. The mean radial velocity of the 
giants, derived from our measurements, is -11.77 km s"', with 
a dispersion of 0.23 km s ' . The mean RV of NGC4349 Nol27 
is -11.40 km s"'; this star is therefore most probably a cluster 
member. 

We use HARPS to observe the NGC 4349 giants since they 
are too faint for Coralie (V = 11.3). To date we have obtained 
about 7 measurements per star spanning ~500 days. The typ- 
ical photon-limited RV precision on these stars amounts to 3- 
4 ms"'. Fig. |2] shows the histogram of the RV dispersion af- 
ter subtraction of the instrumental component. Four stars have 
a moderate scatter of 13-28 m s"', while NGC4349 N0I68 
shows a long-term trend (ctrv - 74 m s"') probably due to a 
stellar companion. Finally, NGC4349 No203 is a short-perio d 
binary exhibiting large RV variations (iMermilhod et al.ll2007l) . 




20 40 60 80 
RV dispersion [m/s] 



20 40 60 80 
RV dispersion [m/s] 



Fig. 2. RV dispersion for giant stars in NGC 2423 and 
NGC 4349. Two stars with dispersions larger than 100 m s ' 
(due to stellar companions) are not shown in these plots. 



NGC4349 No 127 stands out with a scatter of 70 m s^^ and a pe- 
riodic signal, as will be seen in Sect. 4. Using the stable giants, 
we derive an estimated jitter level of 20 m s"' , very similar to the 
NGC 2423 giants. Most interestingly, the high-mass NGC 4349 
stars are therefore suitable for planet search with precise radial 
velocities. 



4. Radial velocity data and orbital solutions 

4.1. NGC2423No3 

We obtained a total of 46 data points for NGC2423 No3 (28 with 
Coralie and 18 with HARPS) over a time span of 1529 days. 
Tables |2] and [3] give the list of these measurements with their 
instrumental error bars. These radial velocities have been com- 
puted using the standard Coralie and HARPS pipelines. As can 
be seen, HARPS measurements are of much higher quality than 
Coralie ones. However, in this case the true uncertainties on the 
stellar radial velocity are dominated by stellar jitter, which we 
have estimated to ~20 m s ' in Sect. 3. In the following we 
therefore quadratically add this value to the instrumental error to 
obtain the final error bar To combine both data sets, we consider 
HARPS and Coralie as two independent instruments, i.e. we in- 
troduce an RV offset between both instruments as a free param- 
eter when fitting a model to the combined data. We note that the 
HARPS-Coralie offset cannot be pre-determined once for good, 
since it depends on the star under consideration, the correlation 
masks used, etc. In the case of NGC2423 No3, the offset is par- 
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B-V [mag] 

Fig. 3. Colour-magnitude diagram for NGC 4349 taken from 
WEBDA. Red giants in our survey are shown as triangles, 
NGC4349 No 127 as a square. 
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Fig. 4. Radial velocity curve as a function of time for 
NGC2423 No3. The fitted orbit corresponds to a planet with 
a minimum mass of 10.6 Mjup and a period of 714 days. Full 
dots indicate Coralie data points, while open dots denote HARPS 
measurements. 



ticularly well constrained since there are time intervals where we 
have data points from both instruments simultaneously. 

A periodic RV variation is clearly visible in the data, with an 
approximate period of 700 days (see Fig.|4]l. Assuming this sig- 
nal is due to an orbiting body (see Sect. 5 for a justification), we 
fit a Keplerian orbit and obtain an orbital period P = 714 days, an 
eccentricity e - 0.21 and a RV semi-amplitude K - 138 m s"', 
leading to a minimum mass mi sin / = 10.6 Mjup and a semi- 
major axis a = 2.10 AU for the companion (see Table |5]l. The 
dispersion of the residuals is 18.3 m s"' and the reduced is 



Table 2. Radial velocity measurements for NGC2423 No3 ob- 
tained with Coralie and their instrumental error bars. All data 
are relative to the solar system barycenter 



JD-2400000 


RV 


Uncertainty 




[kms'l 


[kms'l 


52593.848660 


18.35282 


0.00949 


52683.598863 


18.27881 


0.00849 


52693.581623 


18.26094 


0.01090 


52942.816872 


18.19154 


0.01366 


53016.761041 


18.32692 


0.00918 


53045.537306 


18.32506 


0.01019 


53054.551282 


18.33671 


0.00819 


53098.521695 


18.42907 


0.00835 


53110.534244 


18.42353 


0.01128 


53135.479446 


18.45580 


0.01166 


53140.475994 


18.45805 


0.00814 


53288.866103 


18.38649 


0.01544 


53296.884693 


18.36009 


0.00841 


53330.852860 


18.29186 


0.00849 


53363.767436 


18.23434 


0.00979 


53370.738597 


18.22681 


0.00848 


53444.592134 


18.21294 


0.00762 


53449.544167 


18.20224 


0.01174 


53670.874943 


18.27125 


0.01050 


53700.849066 


18.29006 


0.00711 


53703.852968 


18.29413 


0.00919 


53770.603519 


18.34491 


0.01251 


53776.661243 


18.35691 


0.01013 


54005.871155 


18.36764 


0.01518 


54010.878095 


18.36103 


0.00987 


54029.847346 


18.36718 


0.01191 


54031.843317 


18.33092 


0.00864 


54037.837885 


18.34742 


0.00809 



Table 3. Radial velocity measurements for NGC2423 No3 ob- 
tained with HARPS and their instrumental error bars. All data 
are relative to the solar system barycenter 



JD-2400000 


RV 


Uncertainty 




[km s'] 


[kms'l 


53669.846074 


18.28006 


0.00142 


53674.796790 


18.27805 


0.00147 


53692.862406 


18.27738 


0.00142 


53699.842505 


18.29594 


0.00140 


53721.855948 


18.31043 


0.00127 


53728.752368 


18.31082 


0.00129 


53758.656840 


18.36799 


0.00189 


53764.700939 


18.34264 


0.00107 


53784.639213 


18.38299 
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53817.563693 


18.43727 


0.00142 


53831.568277 


18.41925 


0.00132 


53861.565109 


18.44917 


0.00227 


54050.825005 


18.29213 


0.00148 


54054.864340 


18.31615 
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54078.820633 


18.29799 
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54082.782093 


18.26291 


0.00183 


54114.702923 


18.26573 


0.00296 


54122.715930 


18.25682 


0.00155 



0.84, indicating a good fit to the data, although stellar jitter might 
have been slightly overestimated. 

4.2. NGC4349No127 

We gathered 20 data points for NGC4349 No 127 within a time 
span of 784 days (see Fig.|5]l. All these measurements, listed in 
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Table 4. Radial velocity measurements for NGC4349 No 127 ob- 
tained with HARPS and their instrumental error bars. All data 
are relative to the solar system barycenter. 
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Table |4] were obtained with HARPS and reduced with the stan- 
dard high-precision RV pipeline. As for NGC2423 No3, a jitter 
of 20 m s"' was quadratic ally added to the instrumental error 
bars. Assuming the observed RV variations are due to an orbiting 
body, we fit a Keplerian orbit to the data, which yields an orbital 
period P - 678 days, an eccentricity e = 0.19 and a RV semi- 
amplitude K - 188 m s ' . This corresponds to a minimum mass 
m2 sin i = 19.8 Mjup and a semi-major axis a - 2.38 AU for the 
companion (see Table|5]l. The RMS of the residuals is 12.6 m s ' 
and the reduced 0.52. This unreasonably good value probably 
indicates that stellar jitter was also overestimated in this case. If 
this RV signal is indeed caused by an orbiting body (as suggested 
in Sect. 5), then the companion around NGC4349 No 127 has a 
mass above the D-burning threshold and should therefore be re- 
ferred to as a brown dwarf according to current definitions. Such 
objects are extremely rare around solar-type stars, and it is there- 
fore remarkable to have detected one among a relatively small 
sample of intermediate-mass stars (see discussion in Sect. 6). 

5. Line shape and activity analysis 

It is always necessary to carefully analyze RV measurements of 
giant stars since the RV jitter affecting these stars may induce 
a signal that could be misinterpreted as a planet. The depen- 
dence of RV variability on stellar parameters (mass, metallicity, 
evolutionary stage) is poorly known for giant stars. However, 
there are indications that clump red giants (i.e. in the core He- 
burning phase) are intrinsically more stable than first-ascent 
RGB or AGB stars dBizvaev et al.ll2006h . Short-period (hours to 
days) and long-term (hundreds of days) RV variations ha ve been 
know n to exist in giant sta rs for man y years (see Wal ker et al.l 
ll989i:iLarson et al.il 9931: iHatzes & Cochranlll993lii994 . While 
the short-period modulations are most pr obably due to stellar 
pulsations (e.g. solar- l ike oscillations, see iFrandsen et al.ll2002t 
Ide Ridder et all l2006t iHekker et al.1 l2006h . three mam reasons 
may explain the long-term variations: the presence of an orbit- 
ing body, rotational modulations of surface features and long- 
period, non-radial oscillation modes. To distinguish between 
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Fig. 5. Radial velocity curve as a function of time for 
NGC4349 Nol27. The fitted orbit corresponds to a companion 
with a minimum mass of 19.8 Mjup and a period of 678 days. 



these hypotheses, we use two well-known diagnost ics, the bisec- 
tor velocity span of the cross-correlation function (iHatzesll 19961: 
[Oueloz et al. 2001) and the Call H&K activity index (S-index). 
The bisector span traces line shape variations and should remain 
constant if the measured RVs are due to an orbiting body, while 
the S-index is sensitive to active regions on the stellar surface. 
In the case of rotational modulations, these indicators should ex- 
hibit variations in phase with the radial velocities and the stellar 
rotation period. The rotation period is however difficult to esti- 
mate in our case because relevant parameters such as the stellar 
radius, the projected rotational velocity, the inclination angle and 
possible photometric variability are either not known or too un- 
certain. We therefore rely on the study of bisector and S-index 
measurements to trace inhomogeneities at the stellar surface. 
Finally, long-period, non-radial oscillations in G and K giants 
cannot be excluded but are presently rather hypothetical since 
they have never been unambiguously identified and lack a clear 
theoretical support. If present, such oscillation modes should 
have an effect on the bisector velocity span, although there may 
be cases where this signature is difficult to measure (see for ex- 
ample Hatzes & Cochran 1999, for a detailed discussion). 

Fig. |6] shows the spectrum of NGC2423 No3 and 
NGC4349 No 127 in the Call H region. Although the signal-to- 
noise ratio is low due to the faintness of the stars, no re-emission 
features are visible in the Call H line core. This suggests a low 
level of activity in these giant stars. To compute the Call S-index, 
we closely follow the procedure originally used at Mount Wilson 
(IVaughan et al.|[T978l) . i.e. the fluxes in two narrow bandpasses 
(1 A) centered on the Call H&K cores are integrated and then 
normalized with the integrated flux in two neighbouring spectral 
bands. Given the low signal-to-noise ratio, we carefully checked 
that our measurements and error bars are correct. We made sure 
that the random errors (photon and detector noise) are prop- 
erly propagated and checked that no instrumental effects, such 
as background light pollution, have a significant impact on the 
computed S-index values. The error bars indicate that we are 
able to measure Call flux variations at the 10-20% level. 

Fig. |2] and Fig. [8] show the radial velocity, bisector span and 
S-index for NGC2423 No3 and NGC4349 Nol27 as a function 
of orbital phase. For NGC2423 No3, data from both instruments 
are shown, except the S-index which could not be computed on 



6 C. Lovis and M. Mayor: Planets around evolved intermediate-mass stars 

Table 5. Orbital and physical parameters for the planets/brown dwarfs presented in this paper. 



Parameter 




NGC2423 No3 b 


NGC4349 Nol27 b 


P 


[days] 


714.3 ± 5.3 


677.8 ± 6.2 


T 


[JD-2400000] 


53213 + 21 


54114 ± 34 


e 




0.21 ± 0.07 


0.19 ± 0.07 


V 


[km s-'] 


18.3130 ± 0.0067 


-11.4278 + 0.0118 


(0 


[deg] 


18 + 10 


61 ± 19 


K 


[ms^'l 


137.6 ± 9.1 


188.0 ± 13.0 


a\ sin I 


[10"' AU] 


8.84 


11.5 


f{m) 


[10-«Mo] 


0.180 


0.442 


m-i sin i 


[Mj„p] 


10.6 


19.8 


a 


[AU] 


2.10 


2.38 


N 




46 


20 


Span 


[days] 


1529 


784 


Av (HAIiPS-CoraUe) 


[km s-'] 


-0.0098 ± 0.0067 




a (O-C) 


[ms-'] 


18.3 


12.6 
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Fig. 6. Call H spectral region for NGC2423 No3 (top) and 
NGC4349 Nol27 (bottom). The spectra have been vertically 
shifted for clarity. No re-emission features are visible in these 
giants, indicating a low level of activity. 



Coralie spectra due to the too low signal-to-noise ratio. As can 
be seen for both stars, the bisector span shows no correlation 
with the RV variations. HARPS bisector measurements have dis- 
persions of 3.7 and 12.0 m s"' respectively for NGC2423 No3 
and NGC4349 No 127, which is about one order of magnitude 
smaller than the RV variations. Similarly, the S-index does not 
vary in phase with the RVs (although it shows some variations 
in the case of NGC4349 Nol27). To further study the behaviour 
of the bisector span and the S-index, we computed the Lomb- 
Scargle periodograms of these two quantities, shown in Fig. |9] 
and Fig. [10] To assess the significance of the peaks, we per- 
formed Monte-Carlo simulations to compute their false-alarm 
probabilities (FAPs). No significant periodicities are found, with 
the highest peaks always having at least 50% FAR The 1%-FAP 
level is also indicated in the figures. These results make any fur- 
ther attempts to detect stellar signatures such as rotation or pul- 



Fig. 7. Radial velocity, bisector span and S-index plotted as 
a function of orbital phase for NGC2423 No3. HARPS and 
Coralie measurements are shown as full and open dots, respec- 
tively. No correlations are seen between the RVs and the other 
quantities, supporting the planetary hypothesis. Note the differ- 
ent vertical scales for the RV and bisector span measurements. 



sations in these data very speculative. Finally, we also stress that 
the two stars under consideration are the only ones in their re- 
spective clusters showing such a strong RV signal, despite the 
fact that the other clump giants have similar masses, metallici- 
ties and evolutionary status. Altogether, this leads to the conclu- 
sion that the planetary hypothesis is the best explanation for the 
observed, large-amplitude RV variations. 

6. Discussion and conclusion 

6.1. About planets in open clusters 

Open clusters naturally form homogeneous samples of stars with 
well-constrained basic properties such as mass, metallicity and 
age. As a consequence, searching for planets in open cluster 
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Fig. 8. Radial velocity, bisector span and S-index plotted as a 
function of orbital phase for NGC4349 Nol27. No correlations 
are seen between the RVs and the other quantities, supporting 
the planetary hypothesis. Note the different vertical scales for 
the RV and bisector span measurements. 



Fig. 9. Lomb-Scargle periodograms of the bisector span (top) 
and S-index (bottom) measurements for NGC2423 No3. The 
horizontal dashed line indicates the 1 % false-alarm probability 
threshold. No significant periodicities are found in these quan- 
tities. The frequency of the RV signal is shown as a vertical 
straight line. 



environments has the potential to bring a lot of new informa- 
tion on the statistical properties of planetary systems and con- 
strain planet formation and evolution models. As an example, 
the NGC4349 Nol27 system is extreme under two aspects: it is 
the heaviest star (with an accurate mass determination) around 
which a substellar companion has been detected to date, and it 
is one of the youngest systems known so far. Its age of 200 Myr 
represents an upper limit for the timescale of giant planet/brown 
dwarf formation in the inner regions (a few AUs) surrounding 
intermediate-mass stars. 

With the recent discovery bv lSato et al.l ( l2007l) . we now have 
3 massive planets or brown dwarfs orbiting intermediate-mass 
stars in open clusters. This is most probably only the beginning. 
Our ongoing survey is not complete yet and we plan to extend it 
to other clusters to increase the significance of statistical analy- 
ses. A detailed description of the global results from our survey 
is beyond the scope of this paper and will be presented else- 
where. Incidentally, we also stress the importance of searching 
for planets around FGK dwarfs in open clusters. However, two 
main problems make such a project difficult: the faintness of the 
targets (most clusters have FGK dwarfs fainter than V - 14), 
and stellar activity for clusters younger t han ~1 Gyr (see fo r 
example the RV survey in the Hyades bv lPaulson et al.ll2004l) . 
Nevertheless, future high-precision spectrographs installed on 
8-lOm telescopes and next-generation ELTs should be able to 
carry out RV surveys in a large sample of clusters, thereby 
permitting large-scale comparisons between stars of different 
masses and metallicities. 

There have been a number of studies on the impact 
of cluster environment on planet formation (see for exam- 
pleiA rmitage 2000; Scallv & Clarke 2001; Bonnell et al. 2001; 
ISmiTh & Bonn ell 2001; lMalmberg et al.H2()07l) . Close stellar en- 
counters can potentially destroy, or at least strongly affect, plane- 
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Fig. 10. Lomb-Scargle periodograms of the bisector span (top) 
and S-index (bottom) measurements for NGC4349 No 127. The 
horizontal dashed line indicates the 1 % false-alarm probability 
threshold. No significant periodicities are found in these quan- 
tities. The frequency of the RV signal is shown as a vertical 
straight line. 
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tary systems. Moreover, strong UV radiation from nearby O- and 
B-stars may have a significant impact on the planet formation 
process. Little is presently known about the overall magnitude 
of these effects. Future discoveries of planets in open clusters 
will help put some constraints on these issues. 

6.2. Planet properties as a function of stellar mass 

Although the sample of intermediate-mass stars being searched 
for exoplanets is still limited, we can already try to derive some 
fundamental trends in the characteristics of planetary systems as 
a function of stellar mass. In particular we would like to com- 
pute estimates of the giant planet frequency and typical mass of 
planetary systems for different categories of stellar masses. In 
the following we consider three bins of stellar masses, equally 
spaced on a logarithmic scale: 0.18 < M/Mq < 0.56 (mostly 
M dwarfs), 0.56 < M/Mq < 1.78 (mostly FGK dwarfs) and 
1.78 < M/Mq < 5.62 (intermediate-mass red giants). Precise 
RV surveys targeting the first two bins have been ongoing for 
more than a decade and have yielded reliable esti mates of gi- 
ant planet frequency around sola r-type stars (jJdry et al.l 2007t 
Marcy et al. 2006) and M dwarfs (iBonfils et al.ll2006HEndl et al.l 
2006;: iButler et al.ii2006 ). For the third bin only preliminary re- 



sults are available, but they seem to show an abnormal number 
of massive planets or brown dwarfs compared to the other bins, 
where this kind of objects are extremely rare ("brown dwarf 
desert"). To quantify this more precisely, we have to define a 
region in the parameter space of planetary properties where ob- 
servational biases will not strongly affect the results. In this re- 
spect the limiting bin is obviously the category of intermediate- 
mass stars because of the short duration of the surveys and the 
increased detection limits due to higher stellar masses and RV 
jitter. We therefore choose to focus on massive planets with 
M > 5 Mjup located at orbital distances 0.5 AU < a < 2.5 AU. 
Such planets will always induce RV semi-amplitudes larger than 
45 m s ' on periods shorter than 3 years, and should therefore 
be detectable with ~100% probability by ongoing RV surveys. 
We also set a minimal semi-major axis at a = 0.5 AU because 
intermediate-mass red giants will have engulfed closer-in plan- 
ets during their evolution along the RGB. Obviously, this is only 
a rough estimate since the true minimal semi-major axis at which 
a planet can survive during the RGB phase depends on the de- 
tailed stellar structure and system properties. 

We now have to estimate for each bin how many stars have 
been sufficiently observed by RV surveys to reveal massive 
planets orbiting them. For M dwarfs, we estimate that about 
300 stars meet this criterium if we take into account the high- 
precision surveys c arried out by the HARPS and Elodi e teams 
( Bonfils et al.ll2006) . the C alifornia-Carnegie team ( Butler et al.l 
120061) and the Texas team dEndl et al.ll2"006h . For FGK dwarfs, 
the total number of stars followed by the different groups over 
the past decade amounts to about 3000. For intermedia te-mass 
stars, addin g our survey to the o nes carried out by ISato et al.l 
( l2003h and ISetiawan et al.1 (l2004t) leads us to about 200 tar- 
gets with a mass higher than ~1.8 Mq. Again, these numbers 
are only rough estimates, but they should nevertheless be cor- 
rect enough for a qualitative analysis. Other surveys targeting 
inter mediate-mass stars have recently yielded their first results 
(e.g. lJohnson et al]|2007b and should be included in future anal- 
yses. 

As an input database for exoplanets, we use the Extrasolar 
Planets Encyclopaedia maintained by J. SchneideiQ as of 
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Fig. 11. Average mass of planetary systems as a function of 
stellar mass, taking into account all planets known as of 
February 2007 (see text for details). More massive stars seem 
to harbour more massive planetary systems. 



February 2007. We made two modifications to this database. 
First, we removed the few close binary stellar systems since 
planets in such system s probably have quite different proper- 
ties (Egg enberger et alj 2004). Second, we added to the list a 
few brown dwarf candidates that were not included since their 
minimum masses are significantly higher than 13 Mjup. It is nec- 
essary to take such objects into account since we want an unbi- 
ased census of massive planets and brown dwarfs. However, we 
had to define an upper mass limit to avoid including objects in 
the tail of the mass distribution of stellar companions, since we 
are primarily interested in the distribution of "planetary" com- 
panions, to the left of the brown dwarf desert. We put the limit 
at M = 40 Mjup (in the middle of the desert), which led us to 
add 3 objects to the database (HD 137510b, HD 180777b and 
HD 184860b). We note that the exact value of the cut-off mass 
has no significant impact on the results because of the rarity of 
brown dwarfs. 

We can now compute for each stellar mass bin the fre- 
quency of planets and brown dwarfs in the mass and semi- 
major axis ranges defined above. For M dwarfs, there are no de- 
tected planets meeting the criteria. For solar-type stars, there are 
14 planets out of 3000 stars (0.5%), while for intermediate-mass 
stars 5 planets have been detected among 200 targets (2.5%). 
Assuming the true planet frequency is the same for higher-mass 
stars as for solar-type stars (i.e. 0.5%), there is only a 0.3% prob- 
ability to find 5 planets out of 200 intermediate-mass stars. Even 
if these values are derived from small-number statistics, there 
seems to be a real trend in the sense that more massive stars 
form significantly more massive planets or brown dwarfs than 
lower-mass stars. 

Alternatively, we can also compute for each bin the average 
mass of planetary systems, i.e. the total mass of all planets di- 
vided by the total number of stars in each bin. We choose this 
indicator because it gives useful information on the output of 
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the planet formation process (total mass accreted into planetary 
bodies) and it is relatively free of observational biases. On the 
one hand, it seems reasonable to assume that most of the mass 
in planetary systems is contained in the few largest bodies of 
each system. On the other hand, RV surveys are most sensitive 
to massive planets. We should then be able to obtain reasonable 
estimates for the average planetary system mass in the first two 
bins, while this quantity will probably be underestimated in the 
third bin due to the short duration and the lower sensitivity of the 
surveys. Fig. [TT] shows the results in a histogram, which suggests 
that more massive stars do form more massive planetary systems 
than lower-mass stars, in spite of the previously mentioned ob- 
servational biases. If confirmed, this mass scaling raises ques- 
tions on how to classify objects above 13 Mjup orbiting solar- 
type and intermediate-mass stars. An abrupt transition between 
planets and brown dwarfs would have little meaning if both cat- 
egories of objects are formed by the same physical process. 

Such a scaling in the mass distribution of exoplanets is ex- 
pected in the core-accretion scenario of planet formation since 
more massive stars probably have more massive disks, which 
make it possible to accrete larger amounts of rock, ice and 
gas. However, more quantitative studies are needed. In the disk- 
instability paradigm (e.g. Boss 2006), it is not clear how planet 
formation depends on stellar mass in general, although iBossI 
(1200 6) predicts that this mechanism should not be too sensitive 
to this parameter. It also remains to be seen if high luminosities 
and winds will not prevent the formation of gas giants in the in- 
ner region s surro unding intermediate-mass stars. As an example, 
llda & Liiil ( l2005h predict that the location of the ice boundary at 
larger distances is likely to make the formation process of gas 
giants less efficient. They even predict that the fraction of stars 
harbouring giant planets should decrease beyond 1 Mq, which is 
in contradiction with the results presented in this paper (at least 
for planets heavier than 5 Mjup). The apparently high frequency 
of massive planets around intermediate-mass stars indeed sug- 
gests a rather higher efficiency for the accretion process. 
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